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Abstract 
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Heavy baryons provide us a testing ground to the Standard Model (SM), especially 
to QCD in some aspects. With the accumulation of the experimental data on the heavy 
baryons, more reliable theoretical calculations are needed, although some of them are 
rather complicated. Within the framework of the heavy quark effective theory (HQET) 
which is a model-independent method, the theoretical analysis to the heavy baryons 
containing a single heavy quark is comparatively simple because of the heavy quark 
symmetry [1]. However there are still quantities in this framework which need to be 
determined from nonperturbative QCD. 

QCD sum rule [2], which is regarded as a nonperturbative method rooted in QCD 
itself, has been used successfully to calculate the properties of various hadrons. For 
instances, besides the light mesons [2], light baryons were first considered in Ref. [3]. 
Heavy meson properties were systematically analyzed within the HQET [4,5]. Heavy 
baryons were first discussed in Ref. [6], then masses and Isgur-Wise function for heavy 
baryons were calculated in the HQET to the leading order heavy quark expansion in 
Refs. [7] and [8]. In Ref. [9], the calculation for the heavy baryons began with the 
full theory and results of the calculation were expanded by heavy quark masses. In 
this paper, within the framework of the HQET, we study the heavy baryonic two-point 
correlators to the subleading order of the heavy quark expansion by QCD sum rule 
and obtain results for the heavy baryon masses to that order. 

In the HQET, the heavy quark mass vtiq which is defined perturbatively as the pole 
mass has been removed by the field redefinition. The heavy quark field h v is defined 



by 



P + Q(x) = exp(— irriQV ■ x)h v {x) , 



(1) 



where P + — |(1+ fi). To the order of l/m,Q, the effective Lagrangian for the heavy 
quark is [4] 



£cfr = h v iv ■ Dh v + 



1 



h v {iD) 2 h 



9 T 



h v o uvG^ 1 'h, 



(2) 
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As for the 1/tuq terms, the first one still conserves heavy quark spin symmetry. It is 
the last term which violates the spin symmetry. The heavy baryon mass M is expanded 
as [10] 

5A K 5A G 1 

M = m Q + A + + < sq-1 > +0(— ) , (3) 

m Q m Q m 2 Q 

where A is the heavy baryon mass in the heavy quark limit, which has been calculated 

in Ref. [7]. 5A K and 5A G parameterize the spin-conserved and spin- violated 1/mg 

corrections respectively. All of them characterize the properties of the light degrees 

of freedom, sq denotes the heavy quark spin, and ji stands for the total angular 

momentum of the light degrees of freedom. For Ag baryon, 5A G term vanishes. For 

Eq^ baryons, both ^A^ and 5A G terms are nonvanishing with < sq • ji >= —1 for Eg 

and \ for Eg. 



The heavy baryonic currents j v have been given in Refs. [6] and [7] in the rest 
frame of the heavy baryons. Generally they can be expressed as 

f = e abc (qJ a CTrq b 2 )T'h c v , (4) 

where C is the charge conjugate matrix, r is a flavor matrix, V and V are some gamma 
matrices, and a, b, c denote the color indices. T and V can be chosen covariantly as 

r A = 75 r' A = i , (5) 

for Aq baryon; 

r s = 7" r' s = K + 7^)75 , (6) 

for Eg baryon; 

rv = Y r' s , = -gr + - - + 2 -vW , (7) 

for Eg baryon. The choice of V is not unique. Another kind of baryonic current can 
be obtained by inserting a factor j6 before the T in Eqs. (5-7). The currents given by 
Eqs. (5-7) are denoted as and that with jO insertion as j%. We define the "baryonic 



3 



decay constant" / in the HQET as follows 

<0\f\A Q > = f A u, 

<0\m Q > = fcu, (8) 

<o\j;\x* Q > = ^/ E .u„, 

where u is the spinor and is the Rarita-Schwinger spinor in the HQET respectively. 
fx* is the same as /s in the heavy quark limit. As the mass expansion (3), the square 
of / can be expanded in the same way, 

f 2 = P + J — + J —<s Q - jt > +o(4r) > (9) 

m Q m Q m 2 Q 

where f 2 denotes the leading order result and 5f R2 and 5f° 2 the spin-conserved and 
spin- violated 1/mg corrections respectively. 

The two-point correlator T(u) which we choose for sum rule analyzing in the HQET 



is 



r,» = i J d A xe* kx < 0|T#(aOJ>)|0 >, hi = 1,2, (10) 
where u = 2v ■ k. The hadronic representation of this correlator is 

rv,- (a; = — — -= — H = )a h res. , 11 

where 5 A and <5/ 2 stand for the 1/mg corrections in Eqs. (3) and (9). On the other 
hand, T^u) can be calculated in terms of quark and gluon language with vacuum 
condensates. This establishes the sum rule. We use the commonly adopted quark- 
hadron duality for the resonance part of Eq. (11), 

jm r?f V 



res. 



1 r 00 Im T l (uj') 

IT Juj c Uj' — UJ 



where r? crt (c<j) denotes the perturbative contribution, and u c is the continuum thresh- 
old. In this work, we shall consider only the diagonal correlators (i = j). 

The calculations of T(u) are straightforward. The fixed point gauge is used [11]. 
All the condensates with dimensions lower than 6 are retained. We also include the 
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dimension 6 condensate < q(0)q(x) > 2 in our analysis which is a main contribution. 
We use the gaussian ansatz for the distribution in spacetime for this condensate [12]. 
In the heavy quark limit, we have double checked the analysis of Ref. [7]. We use the 
following values of the condensates, 

< qq > ~ -(0.23 GeV) 3 , 
< a s GG > ~ 0.04 GeV 4 , (13) 
< gcPuvG^q > = ml<qq> , m 2 , ~ 0.8 GeV 2 . 

When u c lies between 2.1 - 2.7 GeV for A Q and between 2.3 - 2.9 GeV for the 
stability window exists. We obtain 

A A = 0.79 ± 0.05 GeV , 

(14) 

fl = (0.3 ±0.1) x 10- 3 GeV 6 , 

for Aq baryon; and 

A s = 0.96 ± 0.05 GeV , 

(15) 

fl = (1.7 ±0.5) x 10- 3 GeV 6 , 

for Sq- 1 baryon. The normalization TrrV = 1 has been used in the analysis. The 
errors quoted in Eqs. (14) and (15) contain only those from the stability of the sum 
rule windows. We still do not know the a s corrections for baryons. Taking the meson 
system as a reference [4], the a s correction is very small for A, however it is very large 
(30%) for /. The numerical results are in agreement with that of Ref. [7], the range 
of the Borel parameter is the same T = 0.4 — 0.7. 

The 1/rriQ corrections to the two-point correlator T(u) can be calculated by in- 
cluding insertions of the 1/wiq operators of the Lagrangian (2) with standard method 
which is shown in Fig. I. The insertions of spin-conserved and spin- violated operators 
are calculated separately. The final form of the sum rules are obtained by performing 
Borel transformation. With some simple tricks [13], the sum rules for the mass and / 
can be separated. The results for the mass of Aq baryon come from the spin-conserved 
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operators only (5A\ = 5A A ): 



T 2 d 

5Aai = -w^ 1 ^' 



T 2 d 
16/1 dT 



(16) 



(5A A2 = -T7T72 3^ /a 2 > 



where 

,2 / ~„ ^2 m 2 



(17) 

and the subscripts 1 and 2 denote j7 and j'2 respectively. The sum rule for / of Aq is 

^ = -5< 1 + 5Er^' (18) 



The masses of baryons Sq and Eg are given in terms of SA K and 5A G . They are 
determined by the following sum rules, 





T 2 d 
1 a K 

16/| dT Sl 


^Af 2 = 


16/| dT Sa 




T 2 d 
16/| dT s 



(19) 



where 



jk = { JL_ r d ^ e -.,T + 3m 2 < m > 2 g 13 < *GG > A/r 
Sl v 2 7 7r 4 5io T 2 5 tt 3 3 ' 

jk _ r d0J0J e e -./T 3m 2 <qq> 2 £ _ 5 < a s GG > 3 2A/T 



jG = < QgGG > T 3 C 2A/T 



(20) 

The sum rules for / are given by 
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It can be seen that while the two diagonal sum rules coincided with each other at the 
leading order, they are no longer the same for the spin-conserved I/ttlq corrections. 

The numerical sum rule results for the I/tuq corrections - 5 A and Sf 2 in Eqs. (3) 
and (9) are given in Tables 1 and 2 and Figs. 2-4. The numerical differences resulting 
from the different choices of j v are not significant. The values of u c are generally smaller 
than the leading order results, but still lie in the allowed range of the leading order 
results. The lower limit of the Borel parameter T = 0.4 is determined by requiring 
that the condensates in Eqs. (17) and (20) have less than 40% contribution. The 
upper limit T = 0.6 is obtained by requiring that the pole contribution is over 70%. 
This window is narrower than the leading order one. In the window T = 0.4 — 0.6 the 
results for 5A A and 5A% are comparatively stable. However from Fig. 4, we see that 
5A^ has no good stability in this window. This is because we have not included the 
Feynman diagrams with internal gluon lines which are expected to be important for 
the spin-violated terms. Therefore the value <5A^ in Table 1 is not reliable. The errors 
quoted in Tables 1 and 2 again only refer to that from the stability of the sum rule 
windows. 

From m Ac and m Ab [14], we determine the heavy quark masses m c = 1.43 ± 0.05 
GeV and m b = 4.83 ± 0.07 GeV. These values give the following results, 

m Sc = 2.52 ± 0.08 GeV , m s * = 2.55 ± 0.08 GeV , (22) 

m S6 = 5.83 ± 0.09 GeV , m s * = 5.84 ± 0.09 GeV . (23) 

From the discussion above we know the individual mass value in Eqs. (22) and (23) 
suffers from the inaccuracy of 5A%. The quantity 

i(m Sq + 2m s * ) = rriQ + A + —(0.22 ± 0.06 GeV 2 ) (24) 

is independent of 5A%, therefore more reliable. It is 2.54 ± 0.08 GeV for c quark case 
and 5.83 ± 0.09 GeV for b quark case. Experimentally m Ec = 2453 ± 0.2 MeV [14]. 
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There is an experimental evidence for S* at m s » = 2530 ± 7 MeV [14]. If we take this 
value for m-£*, we have |(ms c +2ms*) = 2504 ±5 MeV. This is in reasonable agreement 
with the theoretical value. The corresponding quantity for the bottom quark can be 
checked by the experiments in the near future. 

To conclude, we have calculated the 1/mg corrections to the heavy baryon masses 
from the QCD sum rules within the framework of the HQET. This study refines the 
leading order analysis [7]. Furthermore within this framework, we can study the three- 
point correlators which will give the form factors for the weak transitions of the heavy 
baryons [8] to the order of 1/rriQ. It is also viable to include the QCD radiative 
corrections in the leading order and subleading order calculations. Both of these two 
aspects are under our studying. 
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Figure captions 



Fig. 1. The sub leading operator insertions relevant to our analysis. 

Fig. 2. Sum rules for <5A A with (a) j% and (b) ~j v 2 . oo c = 2.0,2.1,2.4 GeV for solid, 
dashed, dash-dotted curves respectively. The sum rule window is T = 0.4 — 0.6 GeV. 

Fig. 3. Sum rules for <SAg with (a) ~j\ and (b) j%. uj c = 2.2, 2.4, 2.7 GeV for solid, 
dashed, dash-dotted curves respectively. The sum rule window is T = 0.4 — 0.6 GeV. 

Fig. 4. Sum rule for <5Ag. The sum rule window is T = 0.4 — 0.6 GeV. 
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Tables 



Table 1. Numerical results for 5 A. 





<5A A (GeV 2 ) 


Mf(GeV 2 ) 


Mg(GeV 2 ) 




cu c =2.1 ±0.1 GeV 


uj c =2A ± 0.2 GeV 






0.09 ±0.03 


0.22 ±0.06 


0.03 ±0.02 


ft 


0.09 ±0.05 


0.21 ±0.06 


0.03 ±0.02 



Table 2. Numerical results for Sf 2 . 





5/ 2 (10- 3 GeV 6 ) 


5/f 2 (10- 3 GeV 6 ) 


5/<? 2 (10- 3 GeV 6 ) 




u c =2.1 ±0.1 GeV 


oo c =2A ± 0.2 GeV 






-0.20 ±0.1 


0.7 ±0.4 


-0.1 ±0.1 


] V 2 


-0.3 ±0.1 


0.8 ±0.5 


-0.1 ±0.1 
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